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§ 1. Introduction 

Photon wave function is a controversial concept. Controversies stem from 
the fact that photon wave functions can not have aU the properties of the 
Schrodinger wave functions of nonrelativistic wave mechanics. Insistence on 
those properties that, owing to pecuharities of photon dynamics, cannot be ren- 
dered, led some physicists to the extreme opinion that the photon wave function 
does not exist. I reject such a fundamentalist point of view in favor of a more 
pragmatic approach. In my view, the photon wave function exists as long as 
it can be precisely defined and made useful. Many authors whose papers are 
quoted in this review share the same opinion and had no reservations about using 
the name the photon wave function when referring to a complex vector-function 
of space coordinates r and time t that adequately describes the quantum state 
of a single photon. 

The notion of the photon wave function is certainly not new, but strangely 
enough it has never been systematically and fully explored. Some textbooks 
on quantum mechanics start the introduction to quantum theory with a dis- 
cussion of photon polarization measurements (cf., for example Dirac [1958][ 



Baym [1969]||Lipkin [IQTSlllCohen-Tannoudji, Diu and Laloe [1977]| ), but in all 



these expositions a complete photon wave function never takes on a specific 
mathematical form. Even Dirac who writes "The essential point is the asso- 
ciation of each of the translational states of the photon with one of the wave 
functions of ordinary wave optics" , never expresses this association in an ex- 
plicit form. In this context he also uses the now famous phrase: "Each photon 
interferes only with itself" which implies the existence of photon wave functions 
whose superposition leads to interference phenomena. 

In the textbook analysis of polarization, only simple prototype two-compo- 
nent wave functions are used to describe various polarization states of the photon 
and with their help the preparation and the measurement of polarization is an- 
alyzed. However, it is not explained, why a wave function should not be used to 
describe also the "translational states of the photon" mentioned by Dirac. After 
such a heuristic introduction to quantum theory, the authors go on to the study 
of massive particles and if they ever return to quantum theory of photons it is al- 
ways within the formalism of second quantization with creation and annihilation 



operators. In some textbooks (cf., for example, Bohm [1954] Power [1964] I one 
may even find statements that completely negate the possibility of introducing 
a wave function for the photon. 

A study of the photon wave function should be preceded by an explanation 
what is the photon and why a description of the photon in terms of a wave func- 
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tion must exist. According to modern quantum field theory, photons, together 
with all other particles (and also quasiparticles, phonons, excitons, plasmons, 
etc.), are the quantum excitations of a field. In the case of photons, these are the 
excitations of the electromagnetic field. The lowest field excitation of a given 
type corresponds to one photon and higher field excitations involve more than 
one photon. This concept of a photon (called the modern photon in a tutorial 
review by Kidd, Ardini and Anton [1989] I enables one to use the photon wave 
function not only to describe quantum states of an excitation of the free field 
but also of the electromagnetic field interacting with a medium. Conceptually, 
the difference between free space and a medium is not essential since the physi- 
cal vacuum is like a polarizable medium. It is filled with all the virtual pairs — 
zero point excitations of charged quantum fields. Therefore, even in free space, 
photons can be also viewed as the excitations of the vacuum made mostly of 
virtual electron-positron pairs fBi alynicki-Birula [1963][ Bjorken [1963] I. 

Even though, in principle, all particles can be treated as field excitations, 
photons are much different from massive particles. They are also different from 
massless neutrinos since the photon number does not obey a conservation law. 
There are problems with the photon localization and as a result the position 
operator for the photon is ill-defined, but the similarities between photons and 
other quantum particles are so ample that the introduction of the photon wave 
function seems to be fully justified and even necessary in order to achieve a 
complete unification of our description of all particles. 



1.1. COORDINATE VS. MOMENTUM REPRESENTATION 

In nonrelativistic quantum mechanics the term coordinate representation is used 
to denote the representation in which the wave function ijj(r) is defined as a 
projection of the state vector on the eigenstates |r) of the components x,y, 
and z of the position operator f , 

m = (rlV-). (1.1) 

The wave function in coordinate representation, therefore, becomes automat- 
ically a function of the eigenvalues of the position operator f . The position 
operators act on the wave function simply through a multiplication. In quan- 
tum mechanics of photons this approach does not work due to difficulties with 
the definition of the photon position operator (cf. IS)). One may still, however, 
introduce functions of the coordinate vector r to describe quantum states of 
the photon. By adopting this less stringent point of view that does not tie the 
wave function in coordinate representation with the formula Hl.l() . one avoids 
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the consequences of the nonexistence of the photon position operator f . In 
principle, any function of r that adequately describes photon states may be 
called a photon wave function in coordinate representation and it is a matter 
of taste and convenience which one to use. It should be pointed out that in 
a relativistic quantum theory, even for particles with nonvanishing rest mass, 
the position operator and the localization associated with it do not live up to 
our nonrelativistic expectations. The differences in localization of photons and, 
say electrons, are more quantitative then qualitative since they amount to the 
"spilling of the wave function" beyond the localization region governed by a 
power law versus an exponential decay. 

The photon wave function in momentum representation has not stirred any 
controversy since the photon momentum operator p is well defined. Its exis- 
tence, as the generator of translations, follows directly from the general theory 



of representations of the Poincare group developed by Wigner [1939] It has 
always been taken for granted by all physicists working in relativistic quantum 
electrodynamics that the notion of the photon wave function in momentum 
representation is well founded. Such wave functions describing initial and fi- 
nal states of photons appear in all formulas for transition amplitudes in the 



S-matrix theory of scattering phenomena (cf., for example, Schweber [1961]| 



Akhiezer and Berestetskii [1965] , B ialynicki-Birula and Bialynicka-Birula [1975]} 
Cohen-Tannoudji, Dupont-Roc and Grynberg [1989]| ). Thus, one may safely as- 



sert that the photon wave function in momentum representation is a well defined 
and fully established object. 

1.2. PHASE REPRESENTATION 

The photon wave functions discussed in this review are distinct from the one- 
mode wave functions that have been introduced in the past (London [1927]} 



Bialynicki-Birula and Bialynicka-Birula [1976] Pegg and Barnett [1988] I to de 



scribe multi-photon states. These functions depend on the phase tp of the field 



and were called the wave functions in the phase representation by Bialynicki-Birula and Bialynicka-Birula [1976] 

The wave functions ^{f) characterize quantum states of a selected mode of the 

quantized electromagnetic field and, in general, they describe a superposition 

of states with different numbers of photons. All spatial characteristics of these 

states are contained in the mode function that defines the selected mode of 

the electromagnetic field. One-mode wave functions vE'(</5) describe properties 

of multi-photon states of the quantized electromagnetic field with all photons 

being in the same quantum mechanical state. In contrast, the photon wave 

function in the coordinate representation can be identified with the mode func- 
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tion itself (cf. tll3|) . It describes a state of a single photon and not a state of 
the quantized field. 



1.3. LANDAU-PEIERLS WAVE FUNCTION 

The concept of the photon wave function in coordinate representation was in- 



troduced for the first time by Landau and Peierls [1930] The same function has 
bee n independently rediscovered more recently by Cook [1982a] Cook [1982b][ 



and Inagaki [1994] The Landau-Peierls proposal has not been met with great 
enthusiasm since their wave function is a highly nonlocal object. 

The nonlocality of the Landau-Peierls wave function is introduced by operat- 
ing on the local electromagnetic field with the integral operator (— A)~^/*, 



((-A)-l/^/)(r) = ^ 



(27r|r-r'|)5/2 



fir'). 



(1.2) 



This integral operator corresponds to a division by \/|k| of the Fourier transform 
and it changes the dimension of the wave function from L~^, characteristic 
of the electromagnetic field, to L"^/^. Therefore the modulus squared of the 
Landau-Peierls wave function has the right dimensionality to be interpreted 
as a probability density to find a photon. In particular, these wave functions 
can be normalized to one with the standard definition of the norm since the 
integral of the modulus squared of the wave function is dimensionless. However, 
as has been already noted by Pauli [1933] despite its right dimensionality the 
nonlocal wave function has serious drawbacks. First, it does not transform under 
Lorentz transformations as a tensor field or any other geometric object. Second, 
a nonlocal wave function taken at a point in one coordinate system depends on 
the values of this wave function in all of space in another coordinate system. 
Third, the probability density defined with the use of a nonlocal wave function 
does not correspond to the probability of interaction of localized charges with 
the electromagnetic field. The vanishing of the wave function at a definite point, 
in Pauli's words (Pauli [1933] I, has "no direct physical significance" because the 
electromagnetic field does act on charges at the points where the probability to 
find a photon vanishes. The Landau-Peierls wave functions can not be used 
as primary objects in the presence of a medium since one is unable to impose 
proper boundary conditions on such nonlocal objects. These functions can be 
introduced, if one wishes so, as secondary objects related to the local wave 
function by a nonlocal transformation {^^. Scalar products and expectation 
values look simpler when they are expressed in terms of the Landau-Peierls 
wave functions but that is, perhaps, their only advantage. 
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1.4. RIEMANN-SILBERSTEIN WAVE FUNCTION 



The mathematical object that fully deserves the name of the photon wave func- 
tion can be traced back to a complexified form of Maxwell's equations that was 
known already at the turn of the century. The earliest reference is the second 
volume of the lecture notes on the differential equations of mathematical physics 
by Riemann that were edited and published by Weber [1901] Various appli- 
cations of this form of Maxwell's equations were given by Silberstein [1907a]j 



Silberstein [1907b]j Silberstein [1914] and Bateman [1915] in the framework of 
classical physics. 

The complex form of Maxwell's equations is obtained by multiplying the first 
pair of these equations 

atD(r,i) = V X H(r,t), V-D(r,i)=0, (1.3) 

by the imaginary unit and then by subtracting from it the second pair 

9tB(r,t) = -V X E(r,i), V-B(r,t)=0. (1.4) 

In the SI units that are used here the vectors D and B have different dimensions 
and prior to subtraction one must equalize the dimensions of both terms. The 
resulting equations in empty space are 



idtFirA) = c:S7 x F(r,i), 
V-F(r,i) = 0, 



where 



F(r,i) = ( 



D(r,i) 



2eo 



2/^0 



(1.5) 
(1.6) 



(1.7) 



an d c = l/^ep/iQ. Around the year 19 30 Majorana (unpublished notes quoted 
by Mignani, Recami and Baldo [1974] \ arrived at the same complex vector ex- 
ploring the analogy between the Dirac equation and the Maxwell equations. 
Kramers [1938] made an extensive use of this vector in his treatment of quantum 
radiation theory. This vector is also a natural object to use in the quaternionic 



formulation of Maxwell's theory ( Silberstein [1914] I. With the advent of spinor 
calculus that superseded the quaternionic calculus, the transformation prop- 
erties of the Riemann-Silberstein vector have become even more transparent. 
When Maxwell's equations were cast into the spinor form ( Laporte and Uhlenbeck [1931] 
Oppenheimer [1931] I, this vector turned into a symmetric second-rank spinor. 
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The use of the Riemann-Silberstein vector as the wave function of the pho- 
ton has been advocated by [Oppenheimer [1931]| [Mohere [1949][ |Good [1957]| 



Bialynicki-Birula [199411 |Sipe [1995]| and 'Bialynicki- Birula [1996a] 

It has aheady been noticed by Silberstcin [1907a] that the important dynam- 
ical quantities associated with the electromagnetic field: the energy density and 
the Poynting vector can be represented as bilinear expressions built from the 
complex vector F. Using modern terminology, one would say that the formulas 
for the energy E, momentum P, angular momentum M, and the moment of en- 
ergy N of the electromagnetic field look like quantum-mechanical expectation 
values 



E = 


I;' 


P = 


2ic 


M = 


2ic 


N = 





(frF*xF, 



(1.8) 
(1.9) 
(1.10) 
(1.11) 



evaluated in a state described by the wave function F. All these quantities are 
invariant under the multiplication of F by a phase factor exp(ia). Such a mul- 
tiplication results in the so called duality rotation (Misner and Wheeler [1957] \ 
of the field vectors 



D'/y^ — cosaD/yeg — sinaB/y^/Io, 
B'/-\/a*o = cosaB/y^ + sinaD/y^. 



(1.12) 
(1.13) 



The coupling of the electromagnetic field with charges fixes the phase a but for 
a free photon one has the same complete freedom in choosing the overall phase 
of the photon wave function as in standard wave mechanics of massive particles. 

The Riemann-Silberstein vector has also many other properties that one 
would associate with a one-photon wave function, except for a somewhat modi- 
fied probabilistic interpretation. Insistence on exactly the same form of the ex- 
pressions for transition probabilities as in nonrelativistic wave mechanics leads 
back to the Landau-Peierls wave function with its highly nonlocal transforma- 
tion properties. 
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§ 2. Wave equation for photons 



The wave equation for the photon is taken to be the complexified form (|1.5|l 
of MaxweU's equations. In order to justify this choice, one may show (cf. 2J 
that the Fourier decomposition of the solutions of this wave equation leads 
to the same photon wave functions in momentum representation that can be 
introduced without any reference to a wave equation directly from the general 
theory of representations of the Poincare group (Bargmann and Wigner |1948|| 
Lomont and Moses [1962 1. There is also a heuristic argument indicating that 
eq! 1)1. 5|l is the risfht choice. Namely, the wave equation 1)1. 5|l can be written 
in the same form as the Weyl equation for the neutrino wave function. As a 
matter of fact, all wave equations for massless particles with arbitrary spin can 
be cast into the same form (cf. tll2ll . 



2.1. WAVE EQUATION FOR PHOTONS IN FREE SPACE 

In order to see a correspondence between Maxwell's equations and quantum me- 



chanical wave equations, one may follow Oppenheimer [1931] and Moliere [1949] 



and rewrite 1)1.5)1 with the use of the spin-1 matrices Sx,Sy, Sz well known from 



quantum mechanics (see, for example, Schiff [1968] I. The matrices that will be 
used here are in a different representation from the one usually used in quan- 
tum mechanics since they act on the Cartesian vector components of the wave 
function and not on the components labeled by the eigenvalues of Sz- That is 
the reason why the matrix Sz is not diagonal. 



(2.1) 



Eq. ()1.5)1 can be written in terms of these spin matrices if the following conver- 
sion rule from vector notation to matrix notation is applied 





" 
















i 




' 


— i 


" 










—i 


, Sy — 













i 













i 





—i 




















a X b = — z(a-s)b. 



The resulting equation 



ihdtF{r,t) = c(s-^V^F{r,t), 



(2.2) 



(2.3) 



is of a Schrodinger type but with a different Hamiltonian. The divergence 
condition 1)1.6)1 can also be expressed in terms of spin matrices either as 



(s.V)s,F = V,F, 



(2.4) 
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or equivalently (Pryce [1948] I as 



V F = AF. 



(2.5) 



The form (|2.3|) of the Maxwell equations compares directly with the Weyl equa- 
tion for neutrinos (Weyl [1929] I 



iWt(j)(r,t) 



c\ cr--V 

i 



(r,t). 



(2.6) 



Eq. (|2.()|l differs from (|2.3|) only in having the Pauli matrices, appropriate 
for spin-1/2 particles, instead of the spin-1 matrices that are appropriate for 
photons. Of course, one may cancel the factors of Ti appearing on both sides 
of eqs. (|2.3|) and (|2.6|) , but their presence makes the connection with quantum 
mechanics more transparent. 

Some authors ( |Oppenheimer [1931][ Ohmura [1956]| Moses [1959] I introduced 
a different, though equivalent, form of the photon wave equation for a four- 
component wave function. The inclusion of the forth component enables one to 
incorporate the divergence condition in a natural way. This approach is directly 
related to the spinorial representation of the photon wave function and will be 
discussed in ^21 

In quantum mechanics the stationary solutions of the wave equation play 
a distinguished role. They are the building blocks from which all solutions 
can be constructed. Stationary solutions of the wave equation are obtained by 
separating the time variable and solving the resulting eigenvalue problem. The 
eigenvalue equation resulting from the photon wave equation H2.3|l is 



c S--V F(r) = fiuj¥{v). 



(2.7) 



Assuming that the photon energy Tiuj is positive, one reads from H2.7(l that the 
projection of the spin on the direction of momentum (helicity) is positive. It 
can easily be checked that one can reverse this sign by changing i into —i in 
the definition H1.7|l of the Riemann-Silberstein vector. Thus, the choice of sign 
in the definition of this vector is equivalent to choosing positive or negative 
helicity, corresponding to left-handed or right-handed circular polarization. In 
order to account for both helicity states of the photons, one has to consider 
both vectors; one may call them F+ and F_. This doubling of the vectors F 



has already been considered by Silberstein [1914] in the context of the classical 
Maxwell equations. 
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Of course, if one is only interested in translating the Maxwell equations into 

a complex form, one can restrict oneself to either F+ or F_. In both cases 
one obtains a one-to-one correspondence between the real field vectors D and 
B and their complex combination F±. However, to have a bona fide photon 
wave function one must be able to superpose diff'erent helicity states without 
changing the sign of the energy (frequency). This can only be done when both 
helicities are described by different components of the same wave function as in 
the theory of spin- 1/2 particles. One can see even more clearly the need to use 
both complex combinations when one deals with the propagation of photons in 
a medium. 



2.2. WAVE EQUATION FOR PHOTONS IN A MEDIUM 

In free space the two vectors F± satisfy two separate wave equations 

i5tF±(r,i) = ±cVxF±(r,i). (2.8) 

In a homogeneous medium, using the values of e and /U for the medium in the 
definition of the vectors F± , 

one also obtains two separate wave equations. The new vectors F± are linear 
combinations of the free-space vectors F j_ , 

= + (2.10) 

F- = ^ [(/f - )n + (i/f + )F°-] . (2.11) 

Thus, the positive and negative helicity states in a medium are certain lin- 
ear superpositions of such states in free space. The necessity to form linear 
superpositions of both helicity states shows up even more forcefully in an inho- 
mogeneous medium, because then it is not possible to split the wave equations 
into two independent sets. For a linear, time- independent, isotropic medium, 
characterized by space-dependent permittivity and permeability, one obtains 
the following coupled set of wave equations 



idtF+{r,t) = v{r){\/ xF+{r,t) 



(2.12) 
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F+(r,i) X Vw(r) 



2v{r) 

idtF_(r,t) = -'y(r)(V X F_(r,t) 
1 



2h{r) 



F_(r,t) X V/i(r)), 



(2.13) 



2,(r)F-(r,t)xV.(r)-^F4r,t)xVMr)), 

where F±(r, t) are built with the values of e(r) and /i(r) in the medium, v{r) = 
l/y'e(r)/i(r) is the value of the speed of light in the medium, and h{r) = 
^/j,(r)/e(r) is the "resistance of the medium" (the sole justification for the use 
of this name is the right dimensionality of Ohm). The divergence condition 
1)1. 6|l in an inhomogeneous medium takes on the form. 



V-F+(r,t): 
V-F_(r,t) 



1 



-F+(r,i)-V^;(r) 



1 



-F_(r,i).V/i(r), 



2v{r) ' ' ^ ' 2h{r) 

^ F_(r,t)-Vw(r) + — ^F+(r,i)-V/i(r). 



2v{r) ' 2h{r) 

The quantum- mechanical form of eqs. (|2.12|) and H2.13|l is 



(2.14) 
(2.15) 



ihdtF+{r,t) = vMr)(s--V)vMr)F+(r,t) 

v/^(s--V)V^F-(i-, t) 
ih^^{^-Vh{v))Y+{v,t). 



(2.16) 



(2.17) 



In view of the coupling in the evolution equations (|2.16(l and H2.17|l between 
the vectors F+ and F_, one has to combine them together to form one wave 
function J-' with six components 



F+ 
F 



(2.18) 



The wave equation for this function can be written in a compact form 



ihdtT{Y,t) = vMr)P3(s--V)V^^-^(r,i) 

I 



(2.19) 
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where the spin matrices Si operate separately on upper and lower components 



S^F+ 



and the three Pauli matrices pi act on !F as follows 



F 

F, 



P2J' = 



-iF_ 
iF, 



P3^ = 



F+ 
F 



(2.20) 



(2.21) 



The divergence conditions (|2.14|l and H2.15|l can also be written in this compact 
form 



V-^(r, t) = ^^(r, t)-Wv{r) + Pii^f(j, t)-Vh{r). (2.22) 

One would not be able to write a linear wave equation for an inhomogeneous 
medium in terms of just one three-dimensional complex vector, without doubling 
the number of components. Note that the speed of light v may vary without 
causing necessarily the mixing of helicities. This happens, for example, in the 
gravitational field (cf. Ullll . It is only the space-dependent resistance h{r) that 
causes mixing. 

It is worth stressing that the study of the propagation of photons in an inho- 
mogeneous medium separates clearly local wave functions from nonlocal ones. 
In free space there are many wave functions satisfying the same set of equations. 
For example, differentiations of the wave functions or integral operations of the 
type H1.2|l do not change the form of these equations. The essential difference 
between various wave functions shows up forcefully in the study of the wave 
equation in an inhomogeneous medium, or in curved space f ijllll . All previous 



studies, except Bialynicki-Birula [1994] were restricted to propagation in free 



space and this very important point was completely missed. The photon wave 
equations in an inhomogeneous medium are not very simple but that is due, 
perhaps, to a phenomenological character of macroscopic electrodynamics. The 
propagation of a photon in a medium is a succession of absorptions and sub- 
sequent emissions of the photon by the charges that form the medium. The 
number of photons of a given helicity is, in general, not conserved in these pro- 
cesses and that accounts for all the complications. The photon wave equations 
in an inhomogeneous medium is describing in actual fact a propagation of some 
collective excitations of the whole system and not just of free photons. 
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2.3. ANALOGY WITH THE DIRAC EQUATION 

The analogy with the relativistic electron theory, mentioned in the Introduction, 
becomes the closest when the photon wave equation is compared with the Dirac 
equation written in the chiral representation of the Dirac matrices. In this 
representation the bispinor is made of two relativistic spinors 



^(r,i) 



(2.23) 



and the Dirac equation may be viewed as two Weyl equations coupled by the 
mass term 



ihdt(j)(r,t) — c(cT- —W)(j>(r,t) + mc'^xi^i't), 
I 

ihdtxi'"^, i) — ^c((T- -V)x(r, t) + TOC^0(r, t). 



(2.24) 
(2.25) 



These equations are analogous to eqs. (|2.1t)|) and H2.17|l for the photon wave 
function. For photons, the role of the mass term is played by the inhomogeneity 
of the medium. 



§3. Photon wave function in coordinate representation 

Despite a formal similarity between the wave equations for the photon H2.16(l and 
H2.17|l and for the electron H2.24(l and H2.25|l . there is an important difference. 
Photons, unlike the electrons, do not have antiparticles and this fact influences 
the form of solutions of the wave equation and their interpretation. 

3.1. PHOTONS HAVE NO ANTIPARTICLES 

Elementary, plane-wave solutions of relativistic wave equations in free space are 
of two types: they have positive or negative frequency, 

exp(— iwt + k-r) or exp(iu;i — k-r). (3-1) 

According to relativistic quantum mechanics, solutions with positive frequency 
correspond to particles and solutions with negative frequency correspond to an- 
tiparticles. Thus, positive and negative frequency parts of the same solution of 
the wave equation describe two different physical entities: particle and antipar- 
ticle. 
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Photons do not have antiparticlcs, or to put it differently, antiphotons are 
identical with photons. Hence, the information carried by the negative frequency 
solutions must be the same as the information already contained in the positive 
frequency solution. Therefore, one may completely disregard the negative fre- 
quency part as redundant. An alternative method is to keep also the negative 
frequency part but to impose an additional condition on the solutions of the 
wave equation. This condition states that the operation of particle-antiparticle 
conjugation 

^'=(r,t)=pi.F*(r,t), (3.2) 

leaves the function J- invariant 

T\v,t)=T{v,t). (3.3) 

This condition is compatible with the evolution equation l|2.19|l and it elimi- 
nates the unwanted degrees of freedom (cf. 'Bialynicki-Birula [1994] ) . Eq. (|3.3|l 
is automatically satisfied if the wave function is constructed according to the 
definition H2.18|l . This follows from the fact that F+ and F_ are complex con- 
jugate to each other. The information carried by the six-component function T 
satisfying the condition H3.3|l is contained in its positive energy part and is the 
same as that carried by the initial Riemann-Silberstein vector F. It follows from 
H3.3|l that the negative frequency part can always be obtained by complex 
conjugation and by an interchange of the upper and lower components of the 
positive frequency part 

.F(-)(r,t) = pi.FW*(r,i). (3.4) 

In this review, I shall use the symbol ^ to denote the properly normalized, 
positive energy (positive frequency) part of the function J- 

*(r,t) =.F(+)(r,t). (3.5) 

This is the true photon wave function. Proper normalization of the photon wave 
function is essential for its probabilistic interpretation and is discussed in SJ3 
Note that the function ^! carries the same amount of information as the original 
Riemann-Silberstein vector since ^ can be constructed from F by splitting this 
vector into positive and negative frequency parts and then using the first part 
as the upper components of and the complex conjugate of the second part as 
the lower components, 



*(r,t) 



■F(+)(r, t)' 
F(-)*(r,t) 



(3.6) 
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The positive frequency part of tlic solutions of wave equations is a well known 
concept also in classical electromagnetic theory where it is called the analytic 



signal ( Born and Wolf [1980] Mandel and Wolf [1995] I . 



3.2. TRANSFORMATION PROPERTIES OF THE PHOTON WAVE FUNC- 
TION IN COORDINATE REPRESENTATION 

In free space, the components of the electromagnetic field form a tensor and 
that allows one to establish the transformation properties of T . Transformation 
properties of the photon wave function are the same as those of T . Under 
rotations, the upper and the lower half of ^' transform as three-dimensional 
vector fields. Under Lorentz transformations, the upper and the lower part also 
transform independently and the corresponding rules can be inferred directly 
from classical electrodynamics (cf., for example, Jackson [1975] \. Under the 
Lorentz transformation characterized by the velocity v, the vectors F-t change 
as follows 



vxF, 7^ v(vF) 

c 7 -I- 1 (? 



(3.7) 



where 7 is the standard relativistic factor 7 — -^/l — j (? . One may check that 
this transformation preserves the square of these vectors. 



(F')2 = (F) 



(3.8) 



This is easily understood if one observes that Fj_ is a combination of the well- 
known scalar invariant S — (eoE^ — B^//io)/2 and the pseudoscalar invariant 
= i/eo/^oE-B of the electromagnetic field 



S±iV. 



(3.9) 



Thus, rotations and Lorentz transformations act on vectors F± as elements of 



the orthogonal group in three dimensions (Kramers [1938] p. 429) leading to 
the following transformation properties of the wave function 



*'(r',i') = 



C 
C* 



*(r,t). 



(3.10) 



where C is a three-dimensional, complex orthogonal matrix and C* is its com- 
plex conjugate. The unification of rotations and Lorentz boosts into one complex 
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orthogonal transformation is even more transparent for infinitesimal transfor- 
mations 

= *+ ((i(57 + p3(5v)-s)^', (3.11) 

where is the vector of an infinitesimal rotation. Under the space reflection 
r — !■ — r, the upper and lower part of ^' do not transform independently but 
are interchanged because D and B transform as a vector and a pseudovector, 
respectively, 

^'i-r,t)^pi^{v,t). (3.12) 

All these transformation properties can also be simply stated in terms of 
second rank spinor fields (cf. H12II . 



3.3. PHOTON HAMILTONIAN 

The operator appearing on the right hand side of the evolution equation H2.19|l 
for the photon wave function is the Hamiltonian operator H for the photon 

H=^)p,{s-^V)^) + n^^p,{s-Vhir)). (3.13) 

In free space this expression reduces to 

Fo = cp3(s--V). (3.14) 

I 

The formulas H3.13|l and (|3.14(l define a Hermitian operator with continuous 
spectrum extending from —oo to oo. Hermiticity is defined here with respect to 
the standard (mathematical) scalar product 

(*i|*2) = y"dV*l(r)*2(r). (3.15) 

Wave functions of physical photons are built from positive-energy solutions 
of the eigenvalue problem for the Hamiltonian, 

J?*(r) = £;*(r). (3.16) 

There is simple relation between positive-energy solutions and negative-energy 
solutions. One may obtain all solutions for the negative energies just by an 
interchange of upper and lower components since piHpi = ~H. Such a simple 
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symmetry of solutions is a result of photons being identical with antiphotons 
and it is not found, in general, for particles whose antiparticles are physically 
distinct. For example, the solutions of the wave functions describing electrons in 
the Coulomb potential of the proton are quite different from the wave function 
of positrons moving in the same potential. In the first case the potential is 
attractive (bound states), while in the second case it is repulsive (only scattering 
states). 

Explicit solutions of the energy eigenvalue problem for photons are easily 
obtained in free space but in the presence of a medium this can be done only in 
special cases. In this respect, wave mechanics of photons is not much different 
from wave mechanics of massive particles, where explicit solutions can also be 
found only for special potentials. 



§ 4. Photon wave function in momentum representation 

The most thorough textbook treatment of quantum mechanics of photons has 



been given by Akhiezer and Berestetskii [1965] who devoted the whole long 



chapter to this problem. Their discussion is limited to momentum represen- 
tation except for a brief subsection under a characteristic title: "Impossibility 
of introducing a photon wave function in the coordinate representation" . This 
impossibility will be addressed in fjSl 

Wave mechanics of photons in momentum representation can be derived di- 
rectly from relativistic quantum kinematics and group representation theory but 
here the analysis will be based on the Fourier representation of the photon wave 
function in coordinate representation. 

In this review I shall use traditionally the wave vector k instead of the photon 
momentum vector p = Kk as the argument of the wave function in momentum 
space. The explicit introduction of Planck's constant is necessary only when 
the proper normalization of the wave function is needed. 



4.1. PHOTON WAVE FUNCTION AS A FOURIER INTEGRAL 

The standard procedure for solving the wave equations H1.5|l or H2.3|l is based on 
Fourier transformation. Since every solution of eqs. H1.5|l and (|l.t)|) is a solution 
of the d'Alembert equation 

(^^-A)F±(r,i)=0, (4.1) 
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the vectors F± can be represented as superpositions of plane waves 

F+{r,t) ^Vhc y^[f+(k)e— *+*-- + r(k)e-*-'''-"-], (4.2) 

F_(r,t) =s/nc /---3 [f_(k)e-*"*+^*'"-+f;(k)e^"*-''*'"-], (4.3) 



where ui = c|k| and the factor ^/fic has been introduced for future convenience. 
The remaining integral has the dimension of 1/length^, so that the Fourier 
coefficients f± have the dimension of length. It has already been taken into 
account in (|4.2(l and (|4.3|l that the vectors F+ and F_ are complex conjugate 
of each other. In order to fulfill Maxwell's equations, the two complex vectors 
f+ (k) and f_ (k) must satisfy the set of linear, algebraic equations that result 
from H1.5|l and H1.6|l . respectively, 

zkxf±(k) = ±|k|F±(k), (4.4) 
k-f±(k) = 0. (4.5) 

Actually, the second equation is superfluous since it follows from the flrst. So- 
lutions of these equations are determined up to a complex factor. Denoting by 
e(k) a normalized solution of the first equation taken with a plus sign 

ikxe(k) = |k|e(k), (4.6) 
e*(k).e(k) - 1, (4.7) 

one can express the vectors f±(k) in the form 

f+(k)=e(k)/(k,l), f_(k) =e*(k)/(k,-l). (4.8) 

The two complex functions /(k. A), where A = ±1, describe the independent 
degrees of freedom of the free electromagnetic field. The vector e(k) can be 
decomposed into two real vectors li(k) that form together with the unit vector 
n(k) — k/|k| an orthonormal set 

e(k) = (li(k) +zl2(k))/N/2, l^(k).l,(k) = (4.9) 
n(k).li(k) = 0, li(k) X l2(k) = n(k). (4.10) 

The only freedom left in the definition of e(k) is its phase. A multiplication 
by a phase factor amounts to a rotation of the vectors li(k) around the vector 
n(k). The same freedom characterizes the coefficient functions /(k. A). This 
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phase may, in general, depend on k and it plays an important role in the study 
of the photon wave function in momentum representation. The final form of the 
Fourier representation for vectors F± is 

F+(r,i) - ^/^y^e(k)[/(k,l)e-^"*+'^•'•'■ + r(k,-l)e'-*-^''•■-], (4.11) 

F_(r,i) =\/^y'^e*(k)[/(k,-l)e-'-*+*-- + r(k,l)e-*-'''--]. (4.12) 

4.2. INTERPRETATION OF FOURIER COEFFICIENTS 

In free space, the energy, momentum, angular momentum, and moment of en- 
ergy of the classical electromagnetic field are given by the expressions p.8|l - 
Hl.ll|l . With the help of the formulas (|4.11l) of (|4.12(l they can be expressed in 



terms of the coefficient functions /(k. A) (cf., for example, Bialynicki-Birula and Bialynicka-Birula [1975] \ 



E = E/(^^^-/*(k.A)/(k,A), (4.13) 
P = E/(^^kr(k,A)/(k,A), (4.14) 



M = Ey (^^/*(k.A)(^^kx-Dk + A?.^)/(k,A), (4.15) 
N = Ey (^^/*(k,A)z;.a;Dk/(k,A), (4.16) 



where 



and 



Dk = 9k + iAa(k), (4.17) 



a(k) = li(k) • akl2(k) - l2(k) • akli(k). (4.18) 
The operation Dk is a natural covariant derivative on the light cone fStaruszkiewicz [1973 



Bialynicki-Birula and Bialynicka-Birula [1975] Bialynicki-Birula and Bialynicka-Birula 



1987] \ 



Note that this operation depends through the vector Q;(k) on the phase conven- 
tion for the polarization vector e(k). It is interesting to note that the vector 
Q;(k) has similar properties to the electromagnetic vector potential. It changes 
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by a gradient under a change of the phase, but its curl is uniquely defined. In- 
deed, it follows from the definition (|4.18() of a(k) and from the orthonormality 
conditions (|4.1Q(I that the vector Q:(k) obeys the equation 

d.,aj - djai = -e^mk/'^^ ■ (4.19) 



This equation determines the Berry phase ( Bialynicki-Birula and Bialynicka-Birula [1987] I 
in the propagation of photons. 

The formulas (|4.13|l and (|4.14|l indicate that /(k, ±1) describe field ampli- 
tudes with energy fiui and momentum ?ik. The formula H4.15|l shows that /(k, 1) 
and /(k, — 1) describe field amplitudes with positive and negative helicity since 
their contribution to the component of angular momentum in the direction of 
momentum is equal to ±?i, respectively. 

The functions /(k, A) have actually a dual interpretation. In classical theory 
they yield full information about the electromagnetic field. In wave mechanics of 
the photon, these functions are the components of the photon wave function in 
momentum representation. In order to distinguish these two cases, I shall denote 
the two components of the photon wave function in momentum representation 
by a new symbol (/)(k, A). Wave function must be normalized and the proper 
normalization of the photon wave function 0(k, A) is discussed in Sj^l 

The expansion of the photon wave function into plane waves has the following 
form 

Sk \ e(k, l)0(k, 1 



(27r)3 



e(k,-l)0(k,-l). 



g-^c.t+^k.r^ (4.20) 



where 



e(k,l)=e(k), e(k, -1) =e*(k). (4.21) 

The integral 14.20|l defines a certain (continuous) superposition of the wave 
functions 0(k, A) with different values of the wave vector. If the photon wave 
function in momentum representation is accepted as a legitimate concept, then 
the superpositions of such functions must also be accepted. Those who are 
not sure about the meaning of a continuous superpositions, in the form of an 
integral, may restrict the electromagnetic field to a box and replace the integral 
by a discrete sum. 



4.3. TRANSFORMATION PROPERTIES OF THE PHOTON WAVE FUNC- 
TION IN MOMENTUM REPRESENTATION 

Transformation rules for the wave function in momentum representation may be 
derived from the transformation properties of the Ricmann-Silberstein vector. 
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They are the same in the classical theory of the electromagnetic field and in 
the quantum theory of photons. Under space and time translations, functions 
(/)(k, A) are multiplied by the phase factors 

0'(k, A) = exp(-iwio + ik-ro)0(k, A), (4.22) 

where (ro,to) is the four-vector of translation. From expressions 1)4.13(1 and 
H4.14|l one may deduce that under rotations and Lorentz transformations, the 
functions (/)(k, A) are also multiplied by phase factors 

0'(k', A) = exp(-iAe(k, A))(j)(k, A), (4.23) 

where the phase function 0(k, A) depends on the Poincare transformation A. 
This transformation property can easily be derived from the transformation law 
of the energy-momentum four-vector for the electromagnetic field. The right 
hand side in the formulas H4.13|l and (|4.14|l has three factors: the integration 
volume d'^A:/(27r)'^|k|, the four-vector (k,uj), and the moduli squared of 0(k, A). 
The integration volume is an invariant (cf., for example, Weinberg [1995][ p. 67) 



and therefore |(/)(k, A)p must also be invariant. An explicit form of the phase 
0(k, A) corresponding to a given Poincare transformation A can be given (cf.. 



for example, Amrein [1969] I but it is not very illuminating. 



§ 5. Probabilistic interpretation 

Probabilistic interpretation of wave mechanics requires, first of all, a definition 
of the scalar product (^'i|5'2) that is to be used in the calculation of transition 
probabilities. The modulus squared of the scalar product of two normalized 
wave functions |(5'i|^'2)P determines the probability of finding a photon in 
the state 'i'l when the photon is in the state ^'2. The probability, of course, 
must be a pure number and — as a true observable — must be invariant under 
all Poincare transformations. The most obvious definition of the scalar product 
H3.15|l can not be used because it is neither Poincare invariant nor dimensionally 
correct. There is essentially only one candidate for the correct scalar product. 
Its heuristic derivation is the easiest in momentum representation. 



5.1. SCALAR PRODUCT 

According to quantum mechanics each photon with momentum Kk carries en- 
ergy hw. Thus, the total number of photons N present in the electromagnetic 
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field is obtained by dividing the integrand in the formula (|4.13|) by hco, 

»=i:/p$hi/*.a)i^ (5^1) 

A photon wave function describes just one photon. The normalized wave func- 
tion must, therefore, satisfy the condition N = 1. Normalized photon wave 
functions in momentum representation satisfy the normalization condition 

The form of the scalar product can be deduced from the expression for the 
norm and it reads 

(*i|*2) =Y.J ^)<^2(k, A). (5.3) 

This scalar product can be also expressed in terms of the photon wave func- 
tions in coordinate representation by inverting the Fourier transformation in 
eq. ijOnjl 



(^(k,A) = ^e*(k,A)- / dVexp(-ik-r)*(r,t), 

\/hc J 



(5.4) 



where the scalar product with e* (k. A) is evaluated for upper and lower com- 
ponents separately and for each A only one of them does not vanish. Upon 
substituting this expression into (|5.3(l . interchanging the order of integrations, 
and using the following properties of the vectors e(k, A) 

e*(k,A).e(k,A') - S^^y, (5.5) 
^e*(k,A)e, (k,A) = (5,, - n,(k)n,(k), (5.6) 

A 

one obtains the following expression for the scalar product in coordinate repre- 
sentation 

The norm associated with this scalar product is 

= ^ /^'"/ d'r'^^'^'^y-hfnr')- (5.8) 
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The scalar product 1)5. 7|l and the associated norm (|5.8|) for photon wave func- 
tions have been arrived at by numerous authors starting from various premises. 



Gross [1964] has proven that this scalar product and this norm are invariant not 



only under Poincare transformations but also under conformal transformations. 



Zeldovich [1965] derived the formula 1)5. 8|l for the number of photons in terms 



of the electromagnetic field vectors. Recently, the norm H5.8f) has been found 



very useful in the formulation of wavelet electrodynamics (Kaiser [1992] ). The 
same expression (|5.H(> for the scalar product can also be derived by consider- 
ing quantum-mechanical expectation values. That approach has been used by 



Good [1957] and is presented below. 



5.2. EXPECTATION VALUES OF PHYSICAL QUANTITIES 

In wave mechanics of photons, the normalized photon wave function (/)(k, A) 
replaces the classical field amplitudes /(k. A). The classical expressions for the 
energy, momentum, angular momentum, and moment of energy become the 
formulas for the quantum-mechanical expectation values 

{E)^{^\H\^), (P> = (vl/jPlvE-), (5.9) 
(M) = (^IJI*), (N) = (*|K|*). (5.10) 

These equations compared with the formulas 1)4.13)1 - 1)4.16)1 enable one to identify 
the operators H, P, J, and K in momentum representation as 

H = huj, (5.11) 
P = hk, (5.12) 

J = kx-Dk + Afi-J^, (5.13) 
I |k| 

K = huj-T>k. (5.14) 

i 

The operators H, P, J, and K are Hermitian with respect to the scalar 
product given by the formula ()5.3)) . 

The formulas 1)5. 11)) ~ ()5. 14)1 are fully consistent with the interpretation of 
0(k, A) as the probability amplitude in momentum representation. The prob- 
ability density to find the photon with the momentum ?ik and the helicity A 
is 

Probability density = -^^-j^j!]-- (5-15) 
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In order to express the quantum-mechanical expectation values (|5.11|I - H5.14|I 
in coordinate representation one must identify the proper form of the scalar 
product for the photon wave function ^. This identification has already been 
made by Good [1957] who compared the classical formula for the energy of 
the electromagnetic field with the quantum-mechanical expression involving the 
Hamiltonian 

{E) - mH\^), (5.16) 

and came to the conclusion that the scalar product for the photon wave function 
has to be modified as follows 

(^il^-a) = ((fr-^\^-^2- (5.17) 
J H 

It is assumed here that the wave functions are built from positive energy states 
only and that guarantees the positive definiteness of the norm 

II^-IP ^ [ dV^^i^-, (5.18) 
J H 

associated with that scalar product. This form of the scalar product leads to 
the following expectation values of the energy, momentum, angular momentum, 
and moment of energy operators 

(E) = Jd^r^'^H-^(^s-^V^^, (5.19) 

(P) = Jd^r^H-^jV^, (5.20) 

(M> ^ Jd^r'i'^H-^r x-V + hs)^, (5.21) 

(N) = [ d^r^'^H'^Hr^. (5.22) 



Thus, the operators H, P, J, and K in coordinate representation have the form 

H = c(s-^v), (5.23) 

P = -V, (5.24) 

i 

J = rx-V + ;is, (5.25) 
i 

K = Hy. (5.26) 
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All these operators preserve the divergence condition and they are Her- 

mitian with respect to the scalar product H5.17|l . It is also reassuring to note 
that the quantum-mechanical operators of momentum and angular momentum 
in coordinate representation have the same form as in standard quantum me- 
chanics. This can be taken as another indication that ^'(r, t) is a legitimate and 
useful object. 

One may prove directly (without using the Fourier expansions) with the help 
of the following identities 

i?P3s/c* = -V^f, iJpsf X s/c^f = (f X -V + ?is)*, (5.27) 

i i 

and with the use of eq. H2.4|l that the expectation values H5.19|l - H5.21|l reduce 
to the classical expressions (ll.8|l " (|l.ll(l when the wave function is replaced by 
the classical electromagnetic field. 

The scalar product H5.17|l in coordinate representation has been obtained 
from the scalar product H5.3|l in momentum representation. However, the scalar 
product that contains the division by the Hamiltonian can be derived on more 
general grounds and its definition does not depend on the choice of representa- 
tion. It has been shown (Segal [1963J Ashtekar and Magnon [1975J I that such 



a scalar product is a general feature of geometric quantization in field theory. 

Even though the number of photons is given by a double integral, so that 
there is no local expression for the photon probability density in coordinate 
space, the expression for the energy has the form of a single integral over |5'(r)p. 
Therefore, one may introduce a tentative notion of the "average photon energy 
in a region of space" and try to associat e a probabilistic interpretation of th e 
photon wave function with this quantity ( Bialynicki-Birula [1994][ Sipe [1995] I. 
More precisely, the quantity psi^), 

, , Ld3r^'t(r)*(r) 

may be interpreted as the probability to find the energy of the photon localized 
in the region In other words, psi^) is the fraction of the average total energy 
of the photon associated with the region il. The probability density p_E(r, t) to 
find the energy of the photon at the point r, 

^t(r,t)vl/(r,t) 
PE(r,t) ^ — , (5.29) 

is properly normalized to one and it also satisfies the continuity equation 

atP£;(r,i)+V-jB(r,t) = 0, (5.30) 



5. PROBABILISTIC INTERPRETATION 



28 



with the normalized average energy flux 

jB(r,i) = ^^^^^^, (5.31) 

as the probability current. The direct connection between the wave function 
\['(r, t) and the average energy density justifies the name "the energy wave func- 



tion" used by Mandel and Wolf [1995]^ It is understandable that the localization 
of photons is associated with their energy because photons do not carry other 
attributes like charge, fermion number, or rest mass. It is worth noting that for 
gravitons not only the probability but even the energy can not be localized (cf., 



for example, Weinberg and Witten [1980] I. The probabilistic interpretation of 



the energy wave function ^ is still subject to all the limitations arising from the 
lack of the photon position operator, as discussed in fJHl In particular, there are 
no projection operators whose expectation values would give the probabilities 

The transition amplitudes, the operators representing important physical 
quantities, and the expectation values can be expressed with equal ease in mo- 
mentum representation and in coordinate representation. For photons moving 
in empty space both representations are completely equivalent and give the same 
results. The only relevant issue is whether a particular superposition of wave 
functions in momentum representation is useful for the description of quantum 
states of the photon. The distinguished and unique feature of the superpo- 
sition given by the Fourier integrals (|4.2()(l is that they represent local fields. 
They have local transformation properties (|3.1U|) and they satisfy local bound- 
ary conditions. Therefore, for photons moving in an inhomogeneous or bounded 
medium, it is the coordinate representation that is preferred because only in this 
representation one may easily take into account the properties of the medium 
(cf. ®. 

5.3. CONNECTION WITH LANDAU-PEIERLS WAVE FUNCTION 

One may easily convert the scalar product (|5.7|l into a standard form containing 
a single integration with the use of the following identity 

1 /■ 3 1 11 
16^7 ''|r-r'|5/2 |r_r"|5/2 ~ |r'-r"|2- ^^'^^^ 

This enables one to convert the double integral H5.7|l into a single integral 

(*i|«'2) = y"dWi(r)$2(r). (5.33) 
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The new functions $ are the Landau-Peierls wave functions and they are related 
to the photon wave functions ^' through the formula 

The form of the scalar product for the Landau-Peierls wave functions is simple 
but one must pay for this simplicity with the nonlocality of the wave func- 
tions. There is also a simple mathematical argument that shows shortcomings 
of the Landau-Peierls wave function. While for every integrable wave func- 
tion 4" the transformation H5.34(l defines the Landau-Peierls wave function <i>, 
the inverse transformation is singular since it contains a nonintegrable kernel 
|r - r'l'^/^ ( |Amrein [1969]! [Mandel and Wolf [1995]] ). This leads to a paradox 



that for many "reasonable" functions <I> (for example, for every function that 
becomes zero abruptly at the boundary) the energy density is infinite. Thus, it 
is much more natural to treat 5* as the primary and $ as the derived object. 



§6. Eigenvalue problems for the photon wave function 

In wave mechanics of photons as in wave mechanics of massive particles, one may 
study eigenvalues and eigenfunctions of various interesting observables. The 
most important observables, of course, are the momentum, angular momentum, 
energy, and moment of energy — the generators of the Poincare group. The 
eigenfunctions of these observables will be given in coordinate representation 
to underscore the validity and usefulness of the photon wave function in this 
representation. 



6.1. EIGENVALUE PROBLEMS FOR MOMENTUM AND ANGULAR MO- 
MENTUM 

The eigenvalue problems for the components of the photon momentum operator 
have the standard quantum-mechanical form 

^'(r) = hh *(r), (6.1) 

and its solutions depend on r through the exponential functions exp(ik-r). 

The eigenvalue problem for the photon angular momentum also has the stan- 
dard quantum-mechanical form. It contains, as usual, the eigenvalue problem 
for the z-component of the total angular momentum 



J2*(r) = M/*(r), 



(6.2) 
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and the eigenvalue problem for the square of the total angular momentum 

J2*(r) = ?iV(J+ l)^'(r). (6.3) 

The solutions of eqs. 1)6.2(1 and H6.3|l are well known vector spherical harmonics 
(cf. , for example, Messiah [1961]| ) . The direct connection between the quantum- 



mechanical eigenvalue problems and multipole expansion in classical electromag- 



netism has been explored systematically for the first time by Moliere [1949] 



6.2. EIGENVALUE PROBLEM FOR THE MOMENT OF ENERGY 

The solution of the eigenvalue problem for the moment of energy shows the ver- 
satility of the calculational methods based on the coordinate representation and 
sheds some light on the problem of the localizability of the photon that is dis- 
cussed in |5| Of course, the same result can be obtained by Fourier transforming 
the solution of the eigenvalue problem obtained in momentum representation. 

The three components of the moment of energy, like the components of angu- 
lar momentum, do not commute among themselves. Therefore, the eigenvalue 
problem can be posed only for one component at a time. Choosing, for definite- 
ness, the z-component, one obtains the following eigenvalue equation 

- ip3(s-V)z* = K*. (6.4) 

The solution of this equation becomes unique when one chooses two additional 
eigenvalue equations to be solved concurrently. For example, eq. (|6.4|l can be 
solved together with the eigenvalue problems for the x and y components of 
momentum since the three operators Px, Py, and Kz commute. The solutions 
for the upper and lower components of the wave function differ only in the sign 
of K and one can solve them independently. When the wave function in the form 
\1/ — ex.p{ikxX + ikyy){'ipx{z),Tpy{z),'ipz{z)) is substituted into (|6.4|l . one obtains 
a set of ordinary differential equations 

ikyZ^i:, - {z^pyY = Kipx, (6.5) 
{zij^x)' -ikxZipz = H^y, (6.6) 
ikxZipy - ikyZipx = Kipz, (6.7) 

where the prime denotes the differentiation with respect to z. These equations 
are solved by the following substitution (/c^ = k'^ + ky) 

d_ 

' dz ' 



i>x = -p2-{kyK -f kx — )ipz, (6.8) 



i d 
V'y = -p2—{-kxK + ky — )%lj^, (6.9) 
k , z dz 
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which resuhs in a Bessel-type equation for ip^, 

z^iP'^ + z^P'^ + {k^ -klz'^)xp,=0. (6.10) 

The physically acceptable solution of this equation is given by the Macdonald 
function of the imaginary index 

/>oo 

^,{z) ^ K,^{k^z) ^ dte-''^"'°'^*cos{Kt). (6.11) 

^0 

The other solution grows exponentially when z — > oo and must be rejected. The 
physical solution falls off exponentially for large \z\ and represents a photon 
state that is localized as much as possible in the z-direction. The remaining two 
components of the eigenfunction are obtained from eqs. (|6.8(l and (|6.9|l . The 
photon wave functions that describe eigenstates of are not normalizable, 
because the spectrum of the eigenvalues is continuous: k can be any real number. 



6.3. PHOTON PROPAGATION ALONG AN OPTICAL FIBER AS A QUAN- 
TUM MECHANICAL BOUND STATE PROBLEM 

The eigenvalue problem for the photon energy operator in the absence of a 
medium is solved by Fourier transformation as described in 21 In ths presence 
of a medium, one can search for eigenstates of the photon Hamiltonian closely 
following the path traveled in nonrelativistic wave mechanics of massive par- 
ticles. This procedure usually involves selecting a set of operators commuting 
with the Hamiltonian and then solving the appropriate set of eigenvalue equa- 
tions. The photon propagation along an infinite cylindrical optical fiber (cf., 



for example, Bialynicki-Birula [1994] I is a good illustration of this approach. In 
order to take care of the boundary conditions at the surface of the fiber, one 
must work in the coordinate representation. 

Consider an infinite, cylindrical optical fiber of diameter a characterized by a 
dielectric permittivity e. The symmetry of the problem suggests the inclusion in 
the set of commuting operators, in addition to the Hamiltonian, the projections 
of the momentum operator and the total angular momentum on the direction 
of the fiber axis. In cylindrical coordinates the eigenvalue equations for the 
z-components of momentum and angular momentum and the Hamiltonian have 
the form 

-ia^* = fc^*, (6.12) 
{-id^ + is-e,))^ = M^-, (6.13) 

-ip3(s-iepdp + -e^d^ + eA))^ ^ (6.14) 
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where Bp, e^, and are the unit vectors along the coordinate hnes and v 
equals to c outside the fiber. Due to the symmetry of the problem, there is 
no coupling between the upper and lower components of and the solution 
of these eigenvalue equations can be sought in the form of a three dimensional 
vector, 



(6.15) 



In order to separate the variables and obtain a set of ordinary differential equa- 
tions, one needs the following differential and algebraic relations 



d^Gp 



(6.16) 
(6.17) 



All unlisted terms of the type (|6.16(l and 16.17|l vanish. The dependence on ip 
and z of all three components tpp, tpip, and ipz of the photon wave function can 
be separated out on the basis of eqs. (|6.12|l and 16.13|l . 

— exp{ikzz) exp{iM(p)f{p). (6.18) 

The three p-dependent components of the wave function satisfy the equations 



M 



'fz + kzfif 



-(*/p), 



-dpfz + kz{ifp) = -/^, 
-dppU-—{ifp) = -/,. 
These equations lead to a Bessel equation for fz 

,2 



d 



1 TTT 
p ' Cip 9 



fz = 0, 



(6.19) 
(6.20) 
(6.21) 

(6.22) 



where fc^ = uP' jv^ — k\. The remaining two functions / can be determined in 
terms of /z. 



/ = ik '^['^^ + kzdp)fz, 
vp ' 



(6.23) 
(6.24) 
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The photon wave function obeys the Bessel equation inside the fiber with one 
value of k± and with a different values of k± in the surrounding free space. 
The behavior of the solution of eq. H6.22|l depends on whether k± is real or 
imaginary. A general solution of this equation is either (for real kj_) a linear 
combination of Bessel functions of the first kind Jm{p) and the second kind 
Ym{p) or (for imaginary k±) a linear combination of modified Bessel functions 
Im{p) and Km{p)- In full analogy with the problem of a potential well in 
quantum mechanics, one can search for bound states in the transverse direction 
by matching a regular oscillatory solution inside (i.e. the Jm{p) function) with 
an exponentially damped solution outside the fiber (i.e. the Km{p) function). 
The matching conditions, well known from classical electromagnetic theory, are 
the continuity conditions for the and field components at the surface of 
the fiber, when p = a. Bound states occur because the speed of light is greater 
in the vacuum than inside the fiber. Therefore, it may happen that k± is real 
inside and imaginary outside the fiber. Since there are two matching conditions 
and only one ratio of the amplitudes inside and outside the fiber, both conditions 
can be satisfied only for a set of discrete eigenvalues of the photon energy Ticu. 
It is worth noting that in order to have an imaginary k± one must have a 
nonvanishing fc^. Thus, a photon may be bound in the plane perpendicular 
to the fiber, but it is always moving freely along the fiber, as in the quantum- 
mechanical description of a charged particle moving in a homogeneous magnetic 
field. This analysis gives an interpretation of electromagnetic evanescent waves 
as quantum bound states. Of course, true bound states of photons, that are 
described by a photon wave function decaying exponentially in all directions, 
are not possible. 



§ 7. Relativistic invariance of photon wave mechanics 

In a relativistically invariant quantum theory, the Poincare transformations are 
represented by unitary operators. The ten Hermitian generators of these trans- 
formations must satisfy the commutation relations characteristic of the Poincare 
group. The ten generators of the Poincare group are identified with the opera- 
tors H, P, J, and K. They generate infinitesimal time translation, space transla- 
tions, rotations, and boosts (special Lorentz transformations), respectively. The 
structure of the Poincare group leads to the following commutation relations 
obeyed by these generators (cf., for example, [Bargmann and Wigner [1948][ 
Bialynicki-Birula and Bialynicka-Birula [1975] Itzykson and Zuber |1980|||Weinberg [1995]| 



8. LOCALIZABILITY OF PHOTONS 



34 



p. 61) 



/ p 




(7.1) 


Jil Jj 


iTiCijkJkt 


(7.2) 


kkj 


— ihcijkKk, 


(7.3) 




= ihc~^6ijH , 


(7.4) 


k,,H 


= ihPi, 


(7.5) 


k,,kj 


= ~ic~'^heijkJk- 


(7.6) 



All the remaining commutators vanish. One may check by a direct calculation 
that the operators H, P, J, and K, given in momentum representation by the 
formulas H5.11|l - H5.14|l and in coordinate representation by the formulas H5.23|l - 
H5.26|l . obey the commutation relation for the generators of the Poincare group. 
In the proof of the commutation relations in momentum representation one 
needs the condition (|4.19|) . Since all generators of the Poincare group are repre- 
sented by operators that are Hermitian with respect to the scalar product (|5.3|l 
or (|5.7I) . the Poincare transformations are represented by unitary operators. 
Therefore, the scalar product is invariant under these transformations and all 
transition probabilities are the same for all observers connected by Poincare 
transformations. Thus, in both coordinate and momentum representations, 
wave mechanics of photons is a fully relativistic theory. 



§8. Localizability of photons 



The problem of localization of relativistic systems has been first posed and solved 
by'Newton and Wigner [1949] and later refined by 'Wight man 1962] According 
to this analysis it is possible to define position operators and localized states for 
massive particles and for massless particles of spin 0, but not for massless parti- 
cles with spin. Thus, the position operator in the sense of Newton and Wigner 
does not exist for photons (cf. also a recent tutorial review on that subject by 
Rosewarne and Sarkar [1992] I. As a simple heuristic explanation, why position 
operator for the photon does not exist, one may observe (Pryce [1948] I that the 
multiplication by r can not be applied to the photon wave function because it 
destroys the divergence condition (|1.6|l . 
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A weaker definition of localization that is applicable even when the position 



operator does not exist, was proposed by Jauch and Piron [1967] and a very 
detailed analysis of this problem has been given by Amrein |1969|| The Jauch- 
Piron localizability allows for noncompatibility of "photon position measure- 
ments" in overlapping regions. The main weakness of such an abstract analysis 
is that an operational definition of the photon position measurement for pho- 
tons has not been incorporated into it. The existence of position measurements 
for photons is just taken for granted regardless of the feasibility of their phys- 
ical realizations. When a realistic model of the photon detector is brought in, 
it is the wave function rather than <I> that appears as the correct probabil- 



ity amplitude for photodetection (Mandel and Wolf [1995]). Thus, in practical 



applications the energy wave functions VE" always seem to play a dominant role. 

It must, however, be stressed that even for massive particles, the localization is 
not perfect, because it is not relativistically invariant. Two observers who are in 
relative motion would not quite agree as to the localization region of a relativistic 
particle. This follows from the fact that the Newton- Wigner wave function 
^NW jg j-giatgfi to the relativistic wave function that transforms locally under 



Poincare transformations by a nonlocal transformation (cf. |Haag [1993] I 



ij^^ir) = J d\'K{r - r')?A(r'), (8.1) 
where the kernel K can be represented in terms of the Macdonald function 

= ]l^i^f^'K,/,imcr/h). (8.2) 

In the limit, when m ^ 0, K{r) 7r/(27rr)~^/^ and (|8.1|l becomes the relation 
(|1.2|l between the local wave function of the photon and the Landau-Peierls wave 
function. Thus, the difference between the localizability of massive particles 
and photons is not that great. In both cases, localization can not be defined 
in a relativistic manner. However, for massive particles departures from strict 
localization are only exponentially small due to the fast decay of the Macdonald 
function in eq. I|8.1|) . In the nonrelativistic limit, when c — *■ oo, the exponential 
tails become infinitely sharp and the localization is restored. 

Difficulties with the position operator for relativistic particles have a pro- 
found origin connected with the structure of the Poincare group. In nonrela- 
tivistic physics the position operator is the generator (up to a factor of mass) 
of Galilean transformations (cf., for example, Gottfried [1966][ Weinberg [1995][ 



p. 62). In a relativistic theory, the Galilean transformations are replaced by the 
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Lorentz transformations and the position operator (multiplied by the mass) is 
replaced by the boost generator K. The main difference between Galilean and 
Lorentz transformation affecting the discussion of localizability is that boost 
generators do not commute. Therefore, one may only hope to localize relativis- 
tic particles in one direction at a time. The possibility to localize photons in 
one direction has been discussed in general terms as the "front" description by 
Acharya and Sudarshan [1960] The eigenfunctions of the boost operator 
given m may serve as an explicit realization of the front description for the 
photon. 

The considerations of photon localizability, while important for the under- 
standing of some fundamental issues, do not influence much the practical ap- 
plications of the photon wave function. All that really should matter there is 
that the wave function be precisely defined and that its interpretation be not 
extended beyond the limits of applicability. 



§ 9. Phase-space description of a photon 

Distribution functions in phase space are a very convenient tool in the descrip- 
tion of statistical properties and the study of the classical limit of wave me- 
chanics. A direct analog of the Wigner function ( jWigner [1932]| ) introduced in 
wave mechanics may also be introduced for photons with the help of the photon 
wave function. This is done by Fourier transforming the product of the wave 
function and its complex conjugate. Fourier transforms of the electromagnetic 
fields similar to the Wigner function have been introduced in optics, first by 



Walther [1968] in the two-dimensional context of radiative transfer theory and 



then by Wolf [1976] and by Sudarshan [1979, 1981a, b] in the three-dimensional 
case. In these papers phase-space distribution functions were defined for the 
stationary states of the electromagnetic field only and they were treated as 
functions of the frequency. The time-dependent distribution functions can be 



defined ( Bialynicki-Birula [1994] I with the use of the time-dependent wave func- 
tion. The only formal difference between the standard definition of the Wigner 
function in nonrelativistic wave mechanics of massive particles and the case of 
photons is the presence of vector indices. Thus, the photon distribution func- 
tion in phase space is not a single scalar function but rather a 6 x 6 Hermitian 
matrix defined as follows 

Wab{v, k, t) = jd'^s e-*-^ ^a{v + s/2, t)-^l{Y - s/2, t). (9.1) 
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Similar multi-component distribution functions arise also for a Dirac particle and 



one can use some of the techniques developed by Bialynicki-Birula Gornicki and Rafelski [1991] 
to deal with such functions. 



Every 6x6 Hermitian matrix can be written in the following block form 



Wab^ 



(9.2) 



where all 3 x 3 matrices are Hermitian. This decomposition can also be 
expressed in terms of the p matrices 

Wab = PoW^j + PiWlj + P2Wfj + pzWl . (9.3) 

The vector indices i and j refer to the components within the upper and lower 
parts of the wave function and the matrices p act on these parts as a whole. The 
most general photon distribution function, as seen from this analysis, is quite 
complicated. In general, when the medium induces mixing of the two polariza- 
tion states, all components of the distribution function are needed. However, 
when photons propagate in free space, only a subset of these components is suf- 
ficient to account for the dynamical properties of photon beams. The simplest 
case is that of a given helicity. A more interesting case is that of an unpolar- 
ized photon beam: a mixture of both helicities with equal weights. This state 
has to be described by the distribution function because a mixed state can not 
be treated by pure Maxwell's theory. In all these cases phase-space dynamics 
can be described by a 3 x 3 Hermitian matrix, i.e., by just one scalar function 
and one vector function. To this end, one may introduce the following reduced 
distribution function 

I^y(r,k,t) = y"d3se-^'^-^i^,(r + s/2,i)V';(r-s/2,t), (9.4) 

where i/^i are the upper or the lower components of the original wave function. 
The Hermitian matrix VFj" can be decomposed into a real symmetric tensor and 
a real vector according to the formula 



T^"=< + ^e.jfe«fc- (9.5) 

The tensor corresponds to the symmetric part of W^" and the vector corresponds 
to the antisymmetric part. The factor of c has been separated out in the second 
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term since the vector u is related to the momentum density, while the trace of 
Wij is related to the energy density ( Bialynicki-Birula [1994] \. 



The equations satisfied by the components of the photon distribution function 
in free space can be obtained from Maxwell's equations H1.5|l and ()1.6|l for the 
vector F, 

dtW.j ^c(k+^v) e^aWij-c(k-^v) Wuemij, (9.6) 

\ Z / m \ Z / m 



k+'-v)w,,=0=(k-^V)W.,. (9.7) 



This leads to the following set of coupled evolution equations for the real com- 
ponents Wij and Ui 

dtWij = -ceukkiWkj - cejikhwki - y (V^Uj + V jUi - SijVkUk), (9.8) 
dtUi = -ceijkkjUk - ^(VjWy - ViWjj), (9.9) 

and to the subsidiary conditions 

ceijkkjUk ^^jWij, (9.10) 
ceijk"^ jUk ^ -ikjWij . (9.11) 

The k-dependent terms in the evolution equations describe a uniform rotation 
of the vector Ui and of the tensor Wij around the wave vector k so that these 
terms can be eliminated by "going to a rotating coordinate frame" . 

With the help of the subsidiary conditions one can eliminate the remaining 
components and obtain from the evolution equations ij^ the equations for w = 
J2 Wit and u 

dtw = -c^ViU^ (9.12) 
dtUi = -2ceijkkjUk-ViW. (9.13) 

These evolution equations do form a simple, self-contained set. However, as it is 
always the case with the phase-space distribution functions in wave mechanics, 
not all solutions of these equations are admissible. Only those distribution 
functions are allowed that can be represented in the form (|9.4|) at the initial 
time (with u and w satisfying the subsidiary conditions (|9.1U|) and (|9.11|l 'l. 
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§ 10. Hydrodynamic formulation 



It has been shown by Madelung [1926] that the Schrodinger wave equation can 



be cast into a hydrodynamic form. In this form, the complex wave function 
is replaced by real variables: the probability density p and the velocity u of 
the probability flow. The wave equation is replaced by the hydrodynamic evo- 
lution equations for the variables p and u. In order to reduce the number of 
functions from four to the original two, one has to impose auxiliary conditions 
— the quantization condition — on the velocity field. Later, other wave equa- 
tions in quantum mechanics (Pauli, Dirac, Weyl) were also presented in the 
hydrodynamic form. The wave equation for the photon wave function is not 



an exception in this respect. It can also be written ( Bialynicki-Birula [1996b] I 
as a set of equations for real hydrodynamic-like variables. Since the Riemann- 
Silberstein vector F carries all the information about the photon wave function, 
one may use F to define these variables. They comprise the energy density p 
and the velocity of the energy flow v, 

p(r,i)=F*(r,t).F(r,i), p(r, t)v(r, i) = ^F*(r, t) x F(r, t), (10.1) 

the components of the following tensor 

U,{v,t) ^ F:{v,t)F,{v,t) + F*{v,t)F,{v,t), (10.2) 

and another vector u, 

p(r,i)u(r,i) = ^(F*(r,t)VF(r,i) - (VF*(r, t))F(r, t)) . (10.3) 

Owing to the existence of the following identities satisfied by the hydrodynamic 
variables 

= 2c, v^Uk = 0, UjUj = - 2i/2, (10.4) 

only one component of tij is arbitrary, but the hydrodynamic equations look 
more symmetric when all the components are treated on equal footing. The 
number of algebraically independent hydrodynamic variables is reduced from 
eight to six (the number of degrees of freedom described by F) by the following 
quantization condition 

j dS- [V X u - -^SiikiviVvj X Vwfe + v^Vtji X Vtu - 2Uiytji x Vvk)] 27rn, 

(10.5) 
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where n is a natural number. This condition must hold for every choice of the 
integration surface and it states, in essence, that the phase of the wave function 
is uniquely defined (up to an overall constant phase). 

The evolution equations for the hydrodynamic variables are 



dtp + {v-V)p 
dtVi + {v-V)vi 



dtUi + {v-V)ui 



-p{V-v), (10.6) 

-dj{~c^p6ij + pviVj + ptij), (10-7) 

1 c 

-{tijVkdkP ~ cSijVkdkP + -z{vidj + v.jdi)p) 
P ^ 

SijVkditki + 2vkdkUj + cEikiUktij + cEjkiUktii 
{vkdktij - UjdkVk) + ^{UkdkVj + tjkdkVi) 

^{vidktkj + Vjdktki + Vkditkj + Vkdjtki), (10.8) 

-^dj [pejki{thmdit„ii + VkdiVi)] . (10.9) 



They must be supplemented by the equations that express the divergence con- 
dition lfO|l 



-^dkiptik) + peijkVjUk 
+ -^{tjkdkUj - Ujdktjk + VkdkVi - VidkVk) = 0, 



4c 
1 



dkipstkm) + ptikUk 



(10.10) 



(10.11) 



2 

+ £; [ejkiiUidkVj - VjdkUi) + SijiitkidkVj - Vjdktki)] = 0. 

Hydrodynamic description of the photon dynamics is not simple but its ex- 
istence underscores again the unification of the quantum theory of the photon 
with the rest of quantum mechanics. Everything that can be done with other 
particles, can be also done for the photons. 



§11. Photon wave function in non-Cartesian coordinate systems 
and in curved space 



It has been shown by Skrotskii [1957] and Plebanski [1960] (for a pedagogical 
review, see Schleich and Scully [1984]| ) that the propagation of the electromag- 
netic field in arbitrary coordinate systems including also the case of curved 
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spacetime may be described by Maxwell's equations, with all the information 
about the spacetime geometry contained in the relations connecting the field 
vectors E, B and D, H. This discovery can be further enhanced by the ob- 



servation ( Bialynicki-Birula [1994] I that in contradistinction to the case of an 
inhomogeneous medium, in the gravitational field the two photon helicities do 
not mix. Xhis follows from tli6 fcict tlia^t for cLrbitrciry metric ^^i/ the constitutive 
relations can be written as a single equation connecting two complex vectors: 
the vector F(r, t) defined by the formula p. 7(1 and a new vector G(r, t) defined 
as 

G(M)^^(5M,M^ („,) 



In curved space (or in curvilinear coordinates) , the constitutive relations for two 
complex vectors F(r, t) and G(r, t) have the form 

= - — {V^g'' + ^goke"'')GJ, (11.2) 
500 

= -^{9^l/V^-^9"'£^kJ)F^, (11.3) 

where g^i, is the metric tensor, g'^'^ is its inverse, and g is the determinant of 
5^1^. The Maxwell equations expressed in terms of vectors G and F in curved 
spacetime are the same as in flat space 

^a^F(r,^) = VxG(r,t), (11.4) 
V-F(r,t) 0, (11.5) 

In these equations, all derivatives are ordinary (not covariant) derivatives as in 
flat space. The whole difference is in the form of the constitutive relations H11.2|l 
and 1)11. 3|l . These relations contain all information about the gravitational field 
or the curvilinear coordinate system. Since the relations between G and F are 
linear, one may write again two separate wave equations for the two helicity 
states as in flat space. By combining these two equations, one obtains the 
following wave equation for the six-component photon wave function .F(r, t) 

idtJ'{r,t)^ P3\7 xg{r,t), (11.6) 

where 

= -^{9v/V^-^P3g°''£^kJ):F\ (11.7) 
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These equations contain only the matrix that does not mix the hchcity states. 

The true photon wave function \l/(r,t) may be introduced only in the time- 
independent case, when the separation of J^(r, t) into positive and negative 
energy parts is well defined. 



§ 12. Photon wave function as a spinor field 

Soon after the formulation of spinor calculus by van der Waerden in the context 



of the Dirac equation, it has been discovered by Laporte and Uhlenbeck [1931] 
that the Maxwell equations can also be cast into a spinorial form. The spinor 
representation of the electromagnetic field and the Riemann-Silberstein vector 
are closely connected. The components of the vector F are related to the com- 
ponents of a second rank symmetric spinor (j)AB 

(t>oo - -F, + iFy, (12.1) 
(j)oi = F,, (12.2) 
cj,n - F,,+tFy, (12.3) 

and the components of the complex conjugate vector F* are related to the 
components of a second rank symmetric primed spinor ^ , 

= -F*-iF;, (12.4) 
= -Ft, (12.5) 
4>^'^' = F*-iF*. (12.6) 

The property that even in curved space both helicities propagate without mixing 
is in the spinorial formalism a simple consequence of the fact that both spinors 
AB and (j)^ ^ satisfy separate wave equations ( Laporte and Uhlenbeck [1931][ 



Penrose and Rindler [1984] ) 



a'^^'^<9^<^Ai3(r,i) = 0, (12.7) 



CA 



,,a^0^'^'(r,i) = 0, (12., 



where the matrices ^ and cr^^,^ are built from the unit matrix and the 
Pauli matrices, 

{a^^'^) = {I,o-)^'^, (12.9) 
Kab') = {I,-^)ab'. (12.10) 
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Under a Lorentz transformation the second rank spinor changes according to 
the formula 



(12.11) 



where Sj^^ is a 2 x 2 matrix of the fundamental (spinor) representation of the 
Lorentz group. Thus, from the point of view of the representation theory of the 
Lorentz and Poincare groups ( |Streater and Wightman [1978]| [Weinberg [1995]| 
p. 231), the photon wave functions for a given helicity are just the three- 
component fields that transform as irreducible representations (1,0) or (0, 1) of 
the proper Lorentz group (without reflections). In order to accommodate reflec- 
tions one must combine both representations and introduce the six-dimensional 
objects T or 

Eqs. 1)12. 7|l (and similarly eqs. p2.8(l ) represent a set of four equations 
obeyed by four components (^n, </>i2, 02i, </'22) of the second order spinor. All 
four equations can be written in the form of the Dirac equation ( Ohmura [1956][ 
Moses [1959]]l 



ihdtMr.t) — a- — V(/)(r, f), 
I 



(12.12) 



where the matrices ot are 
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(12.13) 



One may check that in this formulation the divergence condition takes on the 
following simple algebraic form: (j)i2 = 4>2i- 

The Maxwell equations expressed in spinor notation and the Weyl equation 
provide just the simplest examples from a hierarchy of wave equations for mass- 
less fields described by symmetric spinors 4>BiB-2 - b„ or ^^i-^a ' -^" . All these 
equations have the form ( Penrose and Rindler [1984J I 



CA 



M'/'ABiS2---S„_i(r,t) 

i,A' B'-, B'„ - B. 



0, 

(r,<) = 0. 



(12.14) 
(12.15) 



This universality of massless wave equations for all spins gives an additional 
argument for treating the Riemann-Silberstein vector as the photon wave func- 
tion. 
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§ 13. Photon wave functions and mode expansion of the electro- 
magnetic field 

The concept of the photon wave function is also useful in the process of quan- 
tization of the electromagnetic field. One may simply apply the procedure of 
second quantization to the photon wave function in the same manner as one 
does it for other field operators. In order to see this analogy, one may recall 
that the field operator ?/'(r) for, say the electron field, is built from a complete 
set of wave functions for the electrons (i") ^^'^ from a complete set of wave 
functions for positrons V'n (i") according to the following rule (cf , for example, 
Schweber [1961][ Bialynicki-Birula and Bialynicka-Birula [1975] Weinberg [1995] 



Hr, t) 



^(V^+Cr, t)a„ 



(13.1) 



where a„ and are the annihilation operators for electrons and positrons re- 
spectively. The second part of the field operator is related to the first by the 
operation of charge conjugation performed on the wave functions and on the 
operators. The analog of charge conjugation for photon wave functions is given 
by eq. (|3.2(l . Following this procedure, one may construct the field operator of 
the electromagnetic field in the form 



^(r,i) = ^(*„(r,t)c„ +pi*;(r,t)ct; 



(13.2) 



where the identity of particles and antiparticles for photons has been taken into 
account by using only one set of creation and annihilation operators. The field 
operator (|13.2|) is non-Hermitian but it satisfies the particle-antiparticle conju- 
gation condition (|3.2() . Therefore, it only has six Hermitian components. These 
Hermitian operators are identified as the field operators D(r, t) and B(r, t) and 
they are obtained from J- through the formula 
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(13.3) 



As a direct consequence of the formula (|13.2|l one may identify the photon 
wave functions in the second-quantized theory with the matrix elements of the 
electromagnetic field operators TorTi and B taken between one-particle states 
and the vacuum 



(13.4) 



14. SUMMARY 



45 



In the simplest case of the free field, when the complete set of photon wave 
functions may be labeled by the wave vector k and helicity A, the formula H13.2|l 
takes on the form 



^(r,t) = 



(27r)s 



e(k, 1) (c(k, l)e-"^*+'''-'- + ct(k, - 1) 



jf — zk-r 



e(k, -l)(c(k, -l)e-*'^*+*-'' + ct(k, l)e»'^*-*.r) 



(13.5) 



In the presence of a medium, the expansion H13.2|l of the electromagnetic field op- 
erator requires the knowledge of a complete set of wave functions that satisfy 
the photon wave equation in the medium. These functions are called usually the 

P- 



mode functions of the electromagnetic field (cf., for example, Louisell [1973] 
240 and Mandel and Wolf [19 95]1p. 905) but the term photon wave functions is 



perhaps more appropriate ( jMoses [1973] Bialynicki-Birula and Bialynicka-Birula [1975] 
The advantage of using the terminology of wave functions is that it brings in all 
the associations with wave mechanics and makes the classification of the modes 
more transparent. In particular, one may use the quantum-mechanical notion 
of eigenfunctions and eigenvalues to classify the functions used in the mode ex- 
pansion (Moses [1973] I and also borrow from quantum mechanics the methods 
of proving their completeness ( Bialynicki-Birula and Brojan [1972] I . 

This discussion shows that the photon wave function is not restricted to the 
wave mechanics of photons. The same wave functions also appear as the mode 
functions in the expansion of the electromagnetic field operators. 



§ 14. Summary 

The aim of this review was to collect and explain all basic properties of a cer- 
tain well defined mathematical object — a six-component function of space-time 
variables — that describes the quantum state of the photon. Whether one de- 
cides to call this object the photon wave function in coordinate representation is 
a matter of opinion since some properties known from wave mechanics of mas- 
sive particles are missing. The most essential property that does not hold for 
the photon wave function is that the argument r of the wave function can not 
be directly associated with the position operator of the photon. The position 
operator for the photon simply does not exist. However, one should remem- 
ber that also for massive particles the true position operator exists only in the 
nonrelativistic approximation. The concept of localization associated with the 
Newton- Wigner position operator is not relativistically invariant. Since pho- 
tons can not be described in a nonrelativistic manner, there is no approximate 
position operator. 
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The strongest argument that can be made for the photon wave function in co- 
ordinate representation is based on the most fundamental property of quantum 
states — on the principle of superposition. According to the superposition prin- 
ciple, wave functions form a linear space. By adding wave functions one obtains 
again legitimate wave function. Once this principle is accepted, the existence of 
photon wave functions in coordinate representation follows from the existence 
of the photon wave functions in momentum representation and these functions 
are genuine by all standards; their existence simply follows from relativistic 
quantum kinematics (or more precisely from the representation theory of the 
Poincare group). The Fourier integral (14.20(1 represents a special combination 
of momentum space wave functions with different momenta and as a matter 
of principle, such linear combinations are certainly allowed. One may only ar- 
gue which superpositions to take as more natural or useful but totally rejecting 
the very concept of the photon wave function in coordinate representation is 
tantamount to rejecting the superposition principle altogether. 

There is not much advantage in using the photon wave function in coordinate 
representation to perform calculations for photons moving in free space. The 
relation of this wave function to momentum wave function is so straightforward 
that one may as well stick to momentum representation. It is only in the 
presence of a medium, especially in an inhomogeneous medium, that the photon 
wave function in coordinate representation becomes useful and even essential. 
Only in the coordinate representation one may hope to solve the eigenvalue 
problems and to take into account the boundary conditions. 

The introduction of the wave function for the photon has many significant 
benefits. The photon wave function enables one to formulate a consistent wave 
mechanics of photons that could be often used as a convenient tool in the quan- 
tum description of electromagnetic fields, independently of the formalism of sec- 
ond quantization. In other words, in constructing quantum theories of photons 
one may proceed, as in quantum theory of all other particles, through two 
stages. At the first stage one introduces wave functions and a wave equation 
obeyed by these wave functions. At the second stage one upgrades the wave 
functions to the level of field operators in order to deal more effectively with 
the states involving many indistinguishable particles and to allow for processes 
in which the number of particles is not conserved. Many methods that have 
proven very useful in the study of particles described by the Schrodinger wave 
functions can also be implemented for photons leading to some new insights. 
These methods include relationships between symmetries and operators, the 
definitions of various sets of modes for the electromagnetic field and their com- 
pleteness relations, eigenvalue problems for various observables(|5I), phase-space 
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representation (SQ, and the hydrodynaniic formulation f ijlUI) . Finally, there are 
important logical and pedagogical advantages coming from the use of the photon 
wave function. The quantum mechanical description of all particles, including 
photons, becomes uniform. 
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